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拉丁 超 立 方 抽样 的 自 适 应 高 斯 小 孔 成 像 蝴 蝶 优 化 算法 
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E 要 : ApupHEWEDLAU LIE FAEERE RORROPE E. FRARI, CIO REUREASURUE,. RATAA AUREUS BL 3E SS 
折 小 孔 成 像 蝴 蝶 优 化 算法 。 HEN RU RS NU dM ME 高 种 群 的 多 样 性 ， 

x ; 然后 引入 在 不 同 进化 时 期 自动 调节 搜索 范围 的 自 适 应 最 优 引 导 策 略 ， 平 衡 算 法 的 全 局 和 局 部 搜索 能 

而 提升 算法 的 寻 优 精度 ; 最 后 采用 高 斯 小 孔 成 像 策略 ， 对 最 优 个 体 进行 扰动 ， 使 入 E P A 

一 步 提 升 算法 的 寻 优 精度 并 加 快 算法 的 通过 对 14 个 基准 测试 函数 进行 仿真 实验 以 及 Wilcoxon 秩 和 检验 ， 

结果 表明 改进 算法 的 寻 优 精度 、 收 敛 速度、 稳定 性 和 可 扩展 性 等 性 能 均 得 到 了 较 大 提高 。 

关键 词 : 蝴蝶 优化 算法 ; 拉丁 超 立 方 抽样 ; iE 惯性 权重 ; 高 斯 小 孔 成 像 ; 高 维 优化 
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Self-adaptive Gaussian keyhole imaging butterfly optimization algorithm 
based on Latin hypercube sampling 


Xu Jiet, Lu Haiyan®?t, Zhao Jinjin?, Hou Xinyu?, Lu Mengdie? 
(a. School of Science, b. Wuxi Engineering Technology Research Center for Biological Computing, Jiangnan University, Wuxi 
Jiangsu 214122, China) 


Abstract: Aiming at the shortcomings of Butterfly optimization algorithm, such as poor population diversity, low optimization 
accuracy and slow convergence speed, this paper proposed a self-adaptive Gaussian keyhole imaging butterfly optimization 
algorithm based on Latin hypercube sampling. Firstly, it used a Latin hypercube sampling population initialization strategy to 
enhance the population diversity and thereby improve the overall search ability of the algorithm. Then, it introduced the self- 
adaptive optimal guidance strategy, which can dynamically adjust the search range in different evolutionary periods to balance 
the global and local search capabilities and hence improve the optimization accuracy of the algorithm. Finally, it adopted a 
Gaussian keyhole imaging strategy to disturb the optimal individuals, making the individuals of the population moving close 
to the optimal individuals, so as to further improve the solution accuracy and speed up the convergence of the algorithm. 
Through simulation experiments and Wilcoxon rank sum tests using 14 benchmark functions, the results showed that the 
performance of the improved algorithm is greatly enhanced in terms of optimization accuracy, convergence speed, stability 
and scalability. 

Key words: butterfly optimization algorithm; Latin hypercube sampling; self-adaptive inertia weight; Gaussian keyhole 
imaging; high dimensional optimization 


引言 


优化 算法 一 样 ， 蝴 蝶 优 化 算法 自身 也 存在 着 一 些 不 足 ， 勾 
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分 进化 算法 5 等 ; 内 电 连 接 过 程 算法 四, 模拟 退火 算法 [等 则 是 基于 


中 ,海鸥 优 化 算法 中 等 ， 进化 算法 如 遗传 进化 算法 内 ， 差 


法 的 种 群 多 样 性 差 、 寻 优 精度 不 高 、 收 敛 速度 慢 等 ， 这 坟 
通过 模拟 自然 现象 的 物理 规律 及 自然 界 各 类 生物 种 群 的 。”” 定 程度 上 限制 了 蝴蝶 优化 算法 的 理论 发 展 和 应 用 范围 ,为 此 ， 
活 习 性 和 行为 特点 而 构造 的 一 类 随机 优化 算法 ， 被 称 之 为 ”研究 者 提出 了 不 同 的 改进 策略 。 文 献 [14] 采 用 信息 交叉 共享 机 制 
优化 算法 。 智 能 优化 算法 通常 分 为 三 类 :群体 智能 算法 ， 策略 使 算法 跳 H 
算法 和 基于 自然 现象 的 算法 由。 群体 智能 算法 如 粒子 群 ”新 型 模糊 决策 策略 并 且 构 建 “ 虚 拟 蝴蝶 ”来 提升 算法 的 搜 
力 。 文 献 [16] 通 过 在 花蜜 位 置 引 入 疯狂 因子 以 增加 种 群 多 


局 部 最 优 , 提高 算法 的 收敛 速度 。 文献 [15] 提 


羊 性 ， 获 取 更 好 的 最 优 解 ， 从 而 提高 了 算法 的 寻 优 精度 。 文 


理科 学 等 领域 的 复杂 问题 提供 了 新 的 求解 途径 。 高 了 算法 的 搜索 效率 。 文 献 [18] 通 过 引入 正 余弦 算法 


R 
自然 现象 的 优化 算法 。 这 些 算法 为 解决 大 量 存 在 于 计算 科学 、 工程 。 献 [17] 采 用 动态 调整 转换 概率 策略 以 平衡 全 局 和 局 部 搜索 ， 
= 
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蝴蝶 优化 算法 (butterfly optimization algorithm, BOA)Æ 


[s 


部 算 子 来 提升 种 群 多 样 性 ， 从 而 避免 算法 陷入 局 部 最 优 


化 、 


Brie 


三 


j 于 粒子 滤波 器 优 之 间 的 反馈 共享 交流 网 络 , 提升 了 算法 跳出 局 部 最 优 的 全 
mra 信 感 器 节点 部 轩 、 支持 向 量 机 参数 优化 、 故 ” ”文献 [21] 通 过 引入 混合 惯性 权重 策略 来 提高 种 群 地 多 样 性 
等 实际 应 用 问题 的 研究 中 ! 引 。 但 是 和 其 他 元 启发 式 X 
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年 由 Arora 等 人 9 提出 的 一 种 新 的 智能 优化 算法 。 由 于 ”文献 [19] 将 分 段 权 重 引 入 蝴蝶 算法 中 ， 以 此 来 平衡 全 局 划 


法 因 具 有 原理 简单 、 ， 在 求解 函数 最 优 值 方 男 和 局 部 开发 ， dk 算法 的 寻 优 精度 。 文献 [20] 通 过 构建 贿 
有 较 好 性 能 等 优点 ， 因 此 被 广泛 研究 


EJI o 
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六 献 [22] 采 用 自 适 应 余 切 惯 性 权重 和 逐 维 变 蜡 策略 来 增强 第 
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法 的 自 适 应 性 
SR E s Mg 
来 增强 
上 述 方法 尽管 
量 的 迭代 次 数 ， ERAN UH 
而 且 从 未 将 概率 统计 的 机 


A, F: 拉丁 超 立 方 抽样 的 自 适 应 高 斯 小 了 


需要 


是 升 算 法 的 寻 优 精度 和 收 
阶段 引入 正弦 和 余弦 因子 
的 寻 优 精度 。 
， 但 算法 大 多 都 
E 高 维 函 数 上 人 性 


能 的 相关 研究 较 少 ， 目 关 知识 应 用 到 


寻 优 精度 ， 引入 在 个 同 进化 时 
最 优 引导 策略 ; 最 后 
算法 的 寻 优 精度 并 ] 
出 uic eni 的 1 
度 而 且 进 
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A UL 


羊 的 自 适 应 高 
Y 丁 超 立 方 抽样 种 群 初 
上 局 搜索 能 
， 提 升 算法 的 


斯 小 了 
始 化 策 
力 ; 然后 


mom UE. 


WT FRR Y ESL EL as 


步 提升 


法 的 收敛 速度 。 实 验 表 昌 


本 文 所 提 


的 香味 强 弱 来 调 
蝶 间 的 信息 交互 以 此 来 达到 蝴蝶 


yz 
UuL- 


J 为 有 很 多 不 确定 的 


rH 


， 因 此 Arora uc diis 


建立 了 蝴蝶 优化 算法 : 
度 的 香味 并 且 能 
间 能 够 相互 吸引 。:i 
的 温度 、 papi iR. 


nau 


高 了 算法 的 寻 优 精 


j 蝶 散发 
原 移动 ， 从 而 
立 置 的 更 新 。 
法 在 算法 
化 的 条 件 下 
RUNE REDE ERE 
cui 体 之 


、 自 然 


适应 度 值 决定 。 用 于 4 相关 信息 。 


HRP, H 


1 个 蝴蝶 在 : LAU 


j ERRORS E 
i cin 蝶 位 置 方向 移动 ， 否 贝 
紧 个 体 将 随机 移动 。 
RARUS 随机 移动 的 
索 相 互 转换 的 概率 称 之 为 转 
ERIKU TF: 


优 蝴 蝶 个 体 
翅 嗅 不 到 来 
蝴蝶 个 体 向 最 
阶段 则 称 为 


第 二 ， 假 设 每 只 出 
的 香味 并 朝 记 
自 最 优 蝴蝶 1 
i etant 为 全 局 


1 — 


第 i iE (REDE RU 
的 强度 ，4 是 与 感知 形态 
需 


ARWR RRE 


Q) 


c 是 第 


Ph 刺激 强度 7 是 


AL 


适应 度 最 好 "m 


要 优化 的 函 的 取 值 范围 为 [0 了] 。 
阶段 即 蝴蝶 个 体 是 向 


(ipit pde Dr BER KSN, 


体位 置 移动 表达 式 如 下 


x =d (r xg) fi 


Hr 


N 


Hr 


xX 


芳香 度 ; 


2 


,是 第 i 个 蝴蝶 在 第 1 代 的 位 置 
的 位 置 ， 访 是 第 :个 蝴蝶 在 第 : 代 感 知 到 的 芳香 度 ; 
之 间 的 随机 数 。 

当 r> 忆 时 ,蝴蝶 处 了 
向 是 随机 的 ， 具体 位 


半 中 最 优 蝴 蝶 


r 是 [0,1] 


P. x 是 第 i 个 蝴蝶 在 第 1! 代 的 位 置 
个 蝴蝶 ， 第 上 个 蝴蝶 的 位 
r 是 [0 之 间 的 随机 数 。 


拉丁 超 立 方 抽样 的 自 适 理应 高 斯 小 孔 成 像 蝴 蝶 优 


G3) 


F 局 部 搜索 阶段 即 蝴蝶 个 体 的 移动 方 


(4) 


rh 28 j 
;是 第 i 个 蝴蝶 在 1 代 感 知 到 的 


本 文 对 蝴蝶 算法 进行 


出 了 拉丁 超 立 方 抽样 的 自 


ChinaXiv 合 作 期 刊 


L 成 像 蝴蝶 优化 算法 第 39 卷 第 9 期 


高 斯 小 孔 成 像 蝴蝶 优化 算法 (self-adaptive Gaussian 
keyhole imaging butterfly optimization algorithm based on Latin 
hypercube sampling, SGLBOA)， 该 算法 首先 引入 拉丁 超 立 方 
由 样 初始 化 策略 ， 从 而 提升 算法 的 全 局 搜索 能 力 ， 其 次 ， 引 
入 在 不 同 进化 时 期 自动 调节 搜索 范围 的 自 适应 最 优 引导 策略 


高 算法 的 寻 优 精度 ， 最 后 利用 高 斯 小 孔 成 像 策略 对 最 优 
个 体 进 行 扰动 ， 从 而 进一步 提升 算法 的 寻 优 精度 并 加 快 算法 
的 收敛 速度 。 
拉丁 超 立 方 抽 样 初始 化 策略 


拉丁 超 立 方 抽 样 (Latin hypercube sampling, LHS) 是 
Mckay 55 AU4-T- 1979 年 提出 的 一 种 从 多 元 参数 分 布 中 近似 
随机 抽样 的 方法 。 在 抽样 数 不 多 的 情况 下 ， 随 机 抽样 不 能 很 
好 地 将 样本 分 散 到 整个 区 间 。 与 随机 抽样 不 同 ， 拉 丁 超 立 方 
由 样 具 有 均匀 分 层 的 特性 以 及 可 以 在 较 少 抽样 的 情况 下 ， 得 
尾部 的 样本 值 等 优点 R31。 

算法 的 收敛 速度 和 精度 在 一 定 程度 上 是 受 初始 种 群 好 坏 

影响 的 ， 
时 性 以 及 种 群 的 多 样 性 ， 而 拉丁 超 立 方 抽样 的 均匀 分 层 和 等 
概率 抽样 的 特点 可 以 保证 其 产生 的 变量 覆盖 整个 分 布 空间 ， 
以 将 拉丁 超 立 方 抽样 应 用 于 算法 的 初始 化 可 以 使 初始 种 群 
Bib. (28 sre T8 ZR EE, 以 提高 初始 种 群 的 多 样 性 ， 


随机 产生 的 初始 种 群 无 法 保证 搜索 空间 分 布 的 合 


算法 的 寻 优 性 能 R29。 


ERATA. 
c) 将 变量 X 的 区 间 [bub] 划分 为 N 个 相等 的 子 区 间 。 


于 蝴蝶 优化 算法 的 初始 种 群 是 随机 产生 的 ， 因 此 也 难 
证 搜索 空间 分 布 的 合理 性 以 及 种 群 的 多 样 性 。 为 解决 这 
一 问题 ， 本 文 提出 的 SGLBOA 算法 通过 引入 拉丁 超 立 方 抽 
ii 将 搜索 空 x 间 均 匀 划 分 ， 使 初始 种 群 能 够 较 好 地 禾 盖 整个 
空间 ， 从 而 提高 种 群 的 多 样 性 ， 提 升 算法 的 整体 搜索 能 
w E, 体 步骤 如 下 : 

a) 确 定 种 群 规模 N 和 维 数 D 。 


设 定 变量 x DXF Jy [bub] , 其 中 1b RI ub 分 别 是 变量 的 


d) 在 每 一 维 里 各 个 子 区 间 中 随机 抽取 一 个 点 。 


e) 将 抽取 的 每 一 维 的 点 组 合 形成 初始 种 群 。 

为 了 直观 地 展示 拉丁 超 立 方 抽 样 初始 化 策略 的 有 效 性 ， 
本 文通 过 实验 分 别 采用 随机 初始 化 方式 和 拉丁 超 立 方 抽 样 初始 
化 方式 对 蝴蝶 种 群 进行 初始 化 , 得 到 如 图 1 所 示 的 种 群 分 布 图 。 
图 1(a) 可 以 看 出 ， 采 用 原 BOA 算法 中 随机 方式 初始 
BERT , 种 群 分 布 较为 集中 , 边缘 区 域 蝴蝶 个 体 分 布 较 少 ， 


和 多 样 性 较 差 ， 而 使 用 拉丁 超 立 方 抽样 初始 化 得 到 的 


在 整个 搜索 空间 的 分 布 则 较为 均匀 。 这 从 直观 上 说 明了 
i prn 


初始 化 策略 的 有 效 性 ， 本 文 实验 部 分 也 对 拉 


HK 


样 初始 化 策略 的 有 效 性 进行 了 验证 。 


2.2 自 适 应 最 优 引导 策 略 


” 受 文献 [29] 的 启发 本 文 将 上 述 自 适应 惯性 权重 因子 进 
行 改进 后 引入 蝴蝶 优化 算法 中 以 平衡 算法 的 全 局 和 局 部 搜索 
能 力 ， 改 进 后 的 自 适应 惯性 权重 因子 为 


更 好 地 平衡 鸡 群 算法 的 全 局 和 局 部 搜索 能 力 , 文献 [29] 


引入 了 非 线 性 递减 的 自 适应 惯性 权重 因子 ， 其 表达 式 如 下 : 


In w,4, — IN Whi 
a(t) = t x ——— N Wm 
O - N (5) 


w(t) = en (6) 
w RERA wO 为 第 次 :迭代 的 惯性 权重 因子 ; 
m 为 惯性 权重 因子 的 最 大 值 和 最 小 值 ( we 709 ， 


wE 0) 
In 
a(t)= x crm -Inw,, (8) 
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其 中 : wa =0.9 Wan 704, ERRER H, 权重 因子 较 大 ， 
从 而 使 蝴蝶 个 体 有 较 大 的 搜索 范围 ， 不 仅 能 扩大 搜索 区 域 ， 
还 能 有 利于 跳出 局 部 极 值 ， 而 在 算法 后 期 搜索 阶段 权重 因子 
较 小 ， 有 利于 提高 蝴蝶 群体 的 局 部 搜索 能 力 ， 从 而 有 利于 找 
到 全 局 最 优 值 。 


(a) 随机 初始 


图 1 


化 
两 种 方式 初 


(b) 拉丁 超 立 方 抽样 初始 化 
始 化 种 群 分 布 图 


Fig. 1 Distribution of initialized population in two ways 


蝴蝶 优化 算法 在 全 


局 搜索 式 (4) 中 引入 最 优 蝴蝶 个 体 8* ， 
引导 蝴蝶 向 最 优 个 体 方向 移动 。 受 此 启发 可 知 ， 要 提高 蝴蝶 


优化 算法 的 收敛 速度 和 


同时 为 了 避免 蝴蝶 个 体 陷 入 


w(t) 加 以 平衡 。 从 而 


1+1 
Xi 


对 于 蝴蝶 优化 算法 局 部 


改进 后 的 师 


=w(t)x g" 


随机 性 而 导致 搜索 效率 较 低 和 | 
引入 最 优 蝴蝶 对 


索 阶 段 一 方面 
采用 定向 变异 策略 ， 
n 取 0.9 FF, 
imr", 


ju d Je A 


同时 引入 


前 最 优 值 的 引入 而 导 
r<s7 ， 则 在 局 部 


如 果 


要 敛 精度 ， 需 要 充分 发 挥 最 优 蝴 蝶 个 
Jk g* 的 作用 ， 引 导 种 群 个 体 尽 可 能 地 向 最 优 个 体 的 方向 移动 。 
局 部 最 优 ， 引 入 自 适 应 惯性 权重 


e(t xg ox) 


搜索 阶段 ， 为 避免 算法 因 较 大 的 


蝶 优 化 算法 的 全 局 搜索 公式 为 
(9) 


要 钱 速 度 较 慢 ， 因 此 在 局 部 搜 
明 蝶 个 体 进行 引导 ， 另 一 方面 


判定 系数 7 ， 经 多 次 实验 可 知 
收敛 速度 和 收敛 精度 最 好 ， 故 本 文中 取 w=09 。 
了 搜索 公式 中 引入 扰动 因子 9 以 避免 因 当 
致 算法 在 局 部 搜索 阶段 陷入 
搜索 公式 中 引入 


局 部 极 值 ; 
自 适 应 惯性 权重 因子 


wO 和 黄金 比例 系数 0.618, 使 蝴蝶 个 体 自 适应 的 朝 着 最 优 个 


体 移动 ， 加 快 收敛 速度 。 


当 7>7 时: 


局 部 搜索 阶段 的 具体 表达 式 为 
x =0xg" +(x -xL)x fi (10) 
0=1+ma(1)xtan( rx| 7-2 )) (11) 


MprEn 时: 


xr! 2w(t)x g* +0.618x(r? xx! —xL)x fi 


其 中 ，tmd(D) 是 柯 西 随机 数 ， 


样 性 ; 


r JÆ [0.1] Z^ [R] 


2.3 ”高 斯 小 孔 成 像 策略 


的 随机 数 。 


小 孔 成 像 是 一 种 常见 的 物 ] 


的 一 侧 倾斜 到 另 一 身 


| 上， 在 屏幕 


(12) 
助 于 增加 扰动 因子 2 的 灵活 多 


理光 现象 ， 光 源 通 过 小 孔 从 板 
形成 倒映 像 ， 通 过 调整 小 孔 


屏 与 物体 之 间 的 距离 ， 接 收 屏 上 像 的 大 小 也 会 随 之 发 生化 B31。 


火焰 ?了 在 xX 轴 上 的 投影 为 全 局 最 优 位 置 8* ,坐标 


为 了 使 算法 寻找 到 更 广 的 反 向 位 置 ， 将 小 孔 成 像 原 理 简 
比 到 反 疝 学 习 策 略 中 。 假 设 某 一 空间 中 


个 高 度 为 h 的 


的 上 下 限 


Ha, b (当前 解 空间 的 上 下 限 )。 在 基点 O 上 放置 一 个 有 小 


孔 的 小 孔 ) 


:， 火 焰 通 


x 的 倒映 像 p', 此 时 


8 的 反 向 点 8” ， 如 


过 小 孔 可 以 在 接收 屏 上 得 到 一 个 高 度 为 


图 2 所 示 。 


在 X 轴 得 到 通过 小 孔 成 像 产生 最 优 位 置 


根据 相似 三 角形 原理 可 得 : 


(a+b)/2-8* h 
g” —(a+b)/2 h (13) 


h Hs 
hon, We 的 反 向 解 为 


peth) Ub) e 


2 


(14) 


2n n 


Ch 
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将 惯性 权重 w(?) 以 及 均值 为 0， 方差 为 1 的 高 斯 随机 数 
normrnd(0,1) 9] A Sj e^ 中 ， 得 到 e^ 的 最 


终 反 向 解 为 


(a+b) (a+b) normrmd(0.1)x g* 
TE 


g” =w(t)x 


Eh, n 为 参数 ， 经 多 次 实验 可 知 7 取 2000 时 ,收敛 速度 和 
人 钱 精度 最 好 ， 故 本 文中 取 n=2000。 


图 


2n 


2 小 也 成像 反 向 学 习 示 意图 


(15) 


n 


Fig.2 Schematic diagram of keyhole imaging reverse learning 


对 于 蝴蝶 优化 算法 ， 在 搜索 后 期 
E, 因此 BOA 倾向 于 


集 在 最 优 蝴蝶 的 周 


孔 成 像 反 向 学 习 策 


EE 


WRAAE RR E He R 


略 可 以 增加 种 群 


优化 算法 的 全 局 搜索 能 
2.4 SGLBOA 算法 流程 及 时 间 复 杂 度 分 析 
2.4.1 SGLBOA 算法 流程 


SGLBOA 算法 流程 如 图 


3 所 示 。 


开始 


i 


的 多 样 性 ， 从 而 提升 蝴蝶 


F 早熟 收敛 。 而 引入 小 


设置 


前 法 各 相关 参 歌 , pedro N KERR fax iteratíor 


HERR a, T EE 


刊 用 拉丁 超 立 方 初始 化 第 略 得 到 初始 种 群 


teration ? 


tas GO fos (8) 计 算 自 适应 惯性 积 重 


wD | 


i 


由 式 O) EBEN RO PER 


按 式 (9) 进 行 全 局 位 置 更 新 


由 式 (11) 计 算 抗 动因 子 8 
按 式 (10) 谤 行 局 部 位 置 更 新 


计算 新 的 适应 度 值 , 确定 当前 最 优 值 和 好 


优 位 置 


vy 是 


执行 高 斯 小 孔 威 


像 反 向 学 习 第 略 


根据 式 (2) 更 杂 香 味 浓度 的 感知 形态 


v 


Ap Re CC RHE GERE 


SGLBOA 算法 流程 图 


Flow chart of SGLBOA 


202205.00091 v1 
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2.4.2 SGLBOA 算法 的 时 间 复 杂 度 分 析 


文献 [31] 可 知 蝴蝶 优化 算法 在 种 群 大 小 为 Nw、 维 


D 时 的 时 间 复 杂 度 为 


在 
的 时 间 


度 值 比较 替换 的 时 间 天 、 
ti 解 的 时 间 复 杂 度 下 和 更 新 感觉 


录 最 


T=0(D+/(D) 
SGLBOA 算法 中 , 种 群 大 小 N 、 维 数 D 、 各 参数 设置 
4 、 求 解 适应 度 值 的 时 间 £(D) 、 新 ， 旧 蝴蝶 个 体 适应 


KF 


计算 香味 浓度 的 时 间 复 杂 度 于、 
区 态 阶 段 的 时 间 复 条 度 


记 


T, 都 与 文献 [341 相同。 拉丁 超 立 方 抽样 初始 化 阶段 的 时 间 复 


因此 和 迭代 前 期 的 总 


在 
更 新 公 


只 蝴蝶 进入 全 局 搜索 ， 


T, -O(NxD) 


体 的 时 间 复 杂 度 为 


Tí 2 Ty +0(t +Nf (D)) Max iterx(T, +T; +T) 


=0(D+ f (D) 


局 搜索 阶段 , SGLBOA 算法 在 蝴蝶 个 体 的 全 局 位 置 


式 中 引入 了 自 适应 惯性 权重 wi), 假设: 有 mm(0<sm<N) 


时 间 复 


在 


(0xqXN-m) , 则 满足 
(11) 产 生 扰 动因 子 9 的 时 间 为 ,按照 式 (10) 进 行 位 置 更 新 
度 为 


时 间 为 


= 
更 新 的 


此 
的 时 间 


< 


据 式 (9) 进 行 位 置 的 逐 维 


BR 


式 (7) 和 (8) 产 和 9 


T; 2 O(k +NDk, + Nf (D) ,) - O(D« f (D)) 


局 部 搜索 阶 


段 ， 假 设 满足 ">7 的 蝴蝶 有 4? 只 
"<7 Hyw 


蝶 个 数 为 N-m-g o 


 ， 则 此 时 的 时 间 复 杂 


T; - O(k, + NDks +Nf (D)+t)=0(D+ f (D)) 


满足 "<y 时 ， 假 设 | 
RD. ， 则 此 时 的 时 


E w(7) 的 时 间 为 4 ， 根 
新 的 时 间 为 厂 , 则 全 局 搜索 阶段 的 


假设 


式 
的 


T! -O(NDK,) - O(D) 


外 ， 假 设 
Jj. W 


Tj =0(k,) 
因此 可 得 SGLBOA 算法 的 时 间 复 杂 度 为 


T' =T! + Max, iterx(T/ +T} +T +T} 7) 


-o(D-f(D)) 


于 以 上 分 析 可 知 ,本文 所 提 


算法 相 比 ， 时 间 复 杂 度 并 未 增加 。 
3 ”实验 仿真 及 结果 分 析 


实验 环境 为 windowsl0, 算法 基于 Matlab2020b. 为 了 测 
效 性 ， 引 用 14 个 基准 测试 函 
数 进行 测 斌 分析， 函数 具体 特征 如 表 1 所 示 。 


试 SGLBOA 算法 的 鲁 棒 性 和 


式 (12) 对 蝴蝶 个 体 进行 逐 维 位 置 
|] i Zi BEI 


式 (14) 和 (15) 产 生 £ 的 小 孔 成 像 反 向 解 e" 


的 SGLBOA 算法 和 BOA 


表 1 测试 函数 
Tab. 1 Test functions 
函数 搜索 空间 “特征 ”理论 最 优 值 
Fl: Bentcigar [-10,10] ” 单 峰 0 
F2: Sum of Different Power [100,100] 单 峰 0 
F3: zakharov [-5,10] ” 单 峰 0 
F4: High ConditioneElliptic [-100,100] 单 峰 0 
F5: Sphere [100,100] 单 峰 0 
F6: Matyas [-10,10] ” 单 峰 0 
F7: Griewank [-600,600] 多 峰 0 
F8: ackley [32,2] ”多 峰 0 
F9: Powell [-4,5] 多 峰 0 
F10: Rastrigin [-5.12,5.12] 多 峰 0 
F11: SchaffersF7 [-100,100] 多 峰 0 
F12: Apline [10,10] ”多 峰 0 
F13: Kowalik [-5,5] 4 维 —3.08E-04 
F14: Drop-Wave [-5.12,5.12]. 29€ -1 
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3.1 种 群 规模 N 的 选取 

在 SGLBOA 算法 中 ， 种 群 规模 N 的 不 同 取 值 将 影响 算 
法 的 寻 优 精度 和 收敛 速度 。 为 了 在 算法 中 设置 合适 的 N 值 ， 
选取 单 峰 函 数 Fl 和 F5， 多 峰 函 数 F9 和 F12 对 N 值 进行 敏 
感性 实验 。 实 验 中 分 别 取 WwW=20 、N=30 、N=40 、N=50 和 


D=30 。 实 验 结果 如 图 4 所 示 。 
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(c) F9 
图 4 


(d) F12 
不 同 种 群 规模 函数 收敛 图 
Fig.4 Convergence diagrams with different population sizes 
图 4 可知 ,无 论 是 对 于 单 峰 函数 Fl 和 F5 还 是 对 于 多 
峰 函 数 F9 I F12, 对 于 相同 的 迭代 次 数 ，N=30 时 SGLBOA 
算法 寻 优 精度 最 好 ; 对 于 相同 的 寻 优 精 度 , N=30 时 SGLBOA 
算法 所 需 的 迭代 次 数 最 少 , 综 上 所 述 , 当 种 群 规模 N=30 IST, 
SGLBOA 算法 能 达到 较 好 的 寻 优 精度 和 收敛 速度 ,因此 本 文 
后 续 实验 中 取 N=30 。 
3.2 三 种 改进 策略 的 有 效 性 分 析 

为 了 验证 拉丁 超 立 方 抽样 初始 化 策略 ， 自 适应 最 优 引导 
策略 以 及 高 斯 小 孔 成 像 反 向 学 习 策 略 的 有 效 性 和 可 行 性 ， 本 
文 利 用 上 述 12 个 函数 对 LBOA，LQBOA 以 及 SGLBOA 进 
行 测试 。 这 里 LBOA 是 在 BOA 中 引入 拉丁 超 立 方 抽样 初始 
化 策略 ，LQBOA 是 在 LBOA 中 引入 自 适 应 最 优 引 导 策 略 ， 
SGLBOA 是 在 LQBOA 中 引入 高 斯 小 孔 成 像 反 向 学 习 策略 。 
各 策略 对 应 的 算法 参数 设置 均 与 BOA 相同 ， 即 转换 率 P=0.8 , 窜 
站 数 <= 0.1 ， 初 始 感 知 形态 c =0.01 。 在 W=30 、 Max_iter=100 、 
D-30 的 条 件 下 , 各 算法 均 独 立 运 行 30 次 。 实 验 结果 如 表 2 所 示 。 
表 2 可 知 ，LBOA 在 均值 和 标准 差 方面 整体 上 优 于 
BOA， 这 说 明 加 入 拉丁 超 立 方 抽样 初始 化 策略 提升 了 算法 的 
收敛 精度 和 稳定 性 ， 而 LQBOA 在 均值 和 标准 差 方面 比 
LBOA 更 好 ， 整 体 达 到 了 最 优 ， 这 表明 加 入 自 适 应 最 优 引导 
策略 进一步 提升 了 算法 的 收敛 精度 和 稳定 性 。 对 于 函数 F2、 
F4~F5, F7, F9~F10, F12 和 Fl14, SGLBOA 的 求解 结果 均 为 
理论 最 优 值 0， 同 时 标准 差 也 为 0， 另 外 在 时 间 方 面 ， 
SGLBOA 也 优 于 LBOA 和 LQBOA， 这 表明 高 斯 小 孔 成 像 


策略 进一步 提升 算法 的 寻 优 精度 以 及 算法 的 稳定 性 并 加 
快 算法 的 收敛 速度 。 
3.3 与 其 他 群体 智能 算法 对 比分 析 

为 了 更 全 面 地 验证 SGLBOA 的 有 效 性 ,文中 选取 了 最 近 


两 年 


化 算法 (AOA)B3， 在 函数 寻 优 方面 性 能 较 好 的 鲸鱼 算法 
(WOA)JB4 和 灰 狠 算法 (GWO)B35 
子 群 算法 (PSO) 中 作为 对 比 算法 。 对比 算法 参数 设置 与 其 对 应 1 


R A, S dOTABXLZDAREESQ B iE qu EC dU CERE CT REA 


3r HE LESE GRE REESE Be ET WR S ECS SABHA SEA 


以 及 经 典 


群体 智能 算法 即 粒 


外 ， 对 于 函 


IIT 


所 用 时 间 


| 
^13 | 
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ù 外 ，SGLBOA 均 能 收敛 到 函数 对 应 的 理论 最 优 值 。 对 于 函数 
F8, SGLBOA 也 收敛 到 了 此 函数 常见 的 最 优 值 8.88E-16。 此 
数 F1~F6 fü F8-F10, SGLBOA 收敛 到 j 
相 比 于 其 他 对 比 算法 而 言 者 


里 论 最 优 
LBS. 并 且 除 


函数 F8 


文献 一 致 。 在 N=30 、Max iter=100 、D=30 的 条 件 下 , 各 算法 外 ,均值 和 标准 差 也 都 求解 到 理论 最 优 值 0。 由 此 可 知 
均 独 立 运 行 30 次 。 实 验 结果 如 表 3 所 示 。 SGLBOA 在 寻 优 精度 、 寻 优 速度 方面 以 及 算法 稳定 性 方面 均 
表 3 可 知 ， 对 上 述 测 试 函 数 ， 除 函数 F8、F13 和 Fl14 ”有 较 好 的 效果 。 
表 2 不同 改 进 策略 下 的 实验 结果 
Tab.2 Experimental results of different improving strategies 
m BOA LBOA LQBOA SGLBOA 
Be 均值 标准 差 ”时 间 均值 标准 差 ”时 间 均值 标准 差 ”时 间 均值 标准 差 时间 
Fl 3.32E-05 4.42E-06 0.344  2.52bE-05 2.52E-06 0.385 0.00E+00 0.00E+00 0.473  0.00E-00 . 0.00E-00 0.394 
F2  3.J7E*50  9.38E*50 0.855 5.36E-08 . 3.85E-08 0.8831  0.00E-00 0.00E+00 0.5487  0.00E-00  0.00E-00 0.4616 
F3 7.41E-06 2.32E-06 0.3313  7.36E-06 | 1.27E-06 0.3561  0.00E-00 0.00E+00 0.4178  0.00E-00 ^ 0.00E-00 0.3644 
F4 2.60E-05 3.82E-06 1.0280 | 1.98E-05 2.87E-06 1.0725  0.00E-00 0.00E+00 1.1698  0.00E-00 0.00E+00 0.3342 
F5 1.68E-05 1.80E-06 0.3159 5.13E-06 4.62E-07 0.3418 0.00E+00 0.00E+00 0.4117 0.00E+00 0.00E+00 0.3343 
F6 1.47E-06 7.94E-07 0.2756 2.58E-07 1.67E-07 0.3161 0.00E+00 0.00E+00 0.4172 0.00E+00 0.00E+00 0.3311 
F7 3.69E-04 1.80E-03 0.5548 3.42E-11 5.13E-11 0.6203 0.00E+00 0.00E+00 0.6737 0.00E+00 0.00E+00 0.5910 
F8 1.04E-02 6.30E-03 0.4669 2.92E-04 2.47E-05 0.4440 8.88E-16 8.88E-16 0.4609 8.88E-16 8.88E-16 — 0.3906 
F9 1.37E-05 2.02E-06 0.5732  $8.74E-06 | 1.26E-06 0.6067  0.00E-00 0.00E+00 0.5151  0.00E-00  0.00E-00 0.3920 
F10 1.15E-04 1.15E-04 0.4579 1.06E-05 2.67E-06 0.4353 0.00E+00 0.00E+00 0.4524 0.00E+00 0.00E+00 0.3502 
Fll  5.17E-02 9.10E-03 2.5816 1.54E-02 3.40E-03 2.6321 0.00E+00 0.00E+00 1.7806 0.00E+00 0.00E+00 1.7041 
F12 1.30E-03 2.03E-04 0.3899 2.35E-04 2.07E-05 0.3679  3.71E-176 0.00E+00 0.4553 0.00E+00 0.00E+00 0.3809 
Fl3  1.20E-03 1.20E-03 0.3200 3.96E-04 6.34E-05 0.3452 3.50E-03 3.00E-03 0.4546 6.30E-03 4.20E-03 0.3802 
Fl4 -1.65E+05 5.60E+05 0.5091 -1.02E+00 1.02E-01 0.5713 -1.00E+00 5.19E-06 0.6521 -1.00E+00 $5.19E-06 0.5145 
表 3 与 其 他 群体 智能 算法 的 对 比 实验 结果 
Tab.3 Experimental results and comparison with other swarm intelligence algorithms 
函数 算法 SSA AOA WOA GWO PSO SGLBOA| 函 数 算 法 SSA AOA WOA GWO PSO  SGLBOA 
最 优 值 0.00E+00 4.23E-36 6.73E-10 3.31E+03 2.48E+07 0.00E+00 最 优 值 8.88E-16  8.88E-16 4.75E-12 20.79E+00 20.72E+00 8.88E-16 
均值 2.42E-64  5.78E-09  1.89E-05 1.39E+04 1.24E+08 0.00E+00 均值 8.88E-16 9.36E-04  5.10E-03 20.98E+00 20.99E+00 8.88E-16 
à 标准 差 1.33E-63 1.77E-08 5.78E-05 8.52E+03 4.78E+07 0.00E+00 ü MEX 0.00E-00 — 3.50E-03  2.79E-02  9.28E-02 1.30E-01 0.00E--00 
时 间 /s | 4.4506 0.4493 1.5276 1.3247 0.7644 0.3941 时 间 /s 1.3022 0.4068 0.4429 0.6021 0.6496 0.3906 
最 优 值 0.00E+00 4.87E-202  1.68E-24  6.40E-04 7,19E+17 0.00E+00 最 优 值 0.00E+00 7.18E-36 1.22E-15 31.29E+00  3.25E404 0.00E+00 
均值 2.32E-44 5.02E-16  2.55E-02 3.43E+00 3.46E+23 0.00E+00 均值 3.28E-52  2.30E-03 1.88E-08 6.06E+04 2.06E+05 0.00E+00 
á 标准 差 1.21E-43 2.72E-15  9.28E-02 8.75E+00 1.76E+24 0.00E+00 Á 标准 差 1.90E-51  1.19E-02  9.62E-08  6.23E*04 1.44E+05 0.00E+00 
时 间 /s 1.6254 0.6161 0.7166 0.7319 1.3179 0.4616 时 间 /s 1.3450 0.5121 0.5321 0.6288 0.9692 0.3919 
最 优 值 0.00E+00 2.41E-10 7.59E+04 1.31E*04  2.00E*04 0.00E+00 最 优 值 0.00E+00 0.00E+00 — 1.14E-13 28.62E-00 338.88E+00  0.00E--00 
均值 1.70E-38 450.71E+00 9.82E+04 3.53E+04 3.32E+04 0.00E+00 均值 0.00E+00 5.20E-02  2.21E-07 55.96E+00 540.44E+00 0.00E+00 
m 标准 差 9.29E-38  2.47E*03 1.25E+04 1.04E*04 8.71E+03 0.00E+00 id 示 准 差 0.00E+00 —2.85E-02  9.89E-07 18.70E-00 133.09E-00  0.00E--00 
B]/s 1.2530 0.3500 0.4119 0.4182 0.5587 0.3644  RI/S 1.3079 0.4184 0.4484 0.6006 0.8167 0.3502 
最 优 值 0.00E+00 —3.56E-27 2.49E-11 5.90E-00 7.55E+05 0.00E--00 最 优 值 0.00E-00  L21E-12  1.79E-06 4.33E-01 3.86E+00 0.00E+00 
均值 4.45E-45  7.79E-04  6.69E-06 27.05E-00 3.43E+06 0.00E+00 均值 6.65E-12 4.06E-06  5.96E-05 6.00E-01 5.20E+00 0.00E+00 
i 标准 差 2.44E-44 — 2.90E-03  2.63E-05 15.81E+00 1.60E+06 0.00E--00 S MEZ  3.64E-11 1.48E-05  9.87E-05  1.27E-01 5.61E-01 0.00E--00 
Hs — 1.6806 0.7465 0.8140 0.8764 1.4982 1.0860 时 间 /s ^ 2.6035 1.5506 1.7544 1.6437 2.4154 1.7041 
最 优 值 0.00E+00 1.70E-50 3.38E-15 2.90E-03  54.46E*00 0.00E+00 最 优 值 0.00E+00 9.32E-25 5.97E-11  6.14E-01  70.31E*00 0.00E+00 
均值 1.27E-56  L70E-04 1.51E-11 1.52E-02 144.17E+00 0.00E+00 均值 1.02E-24 2.14E-06 3.02E-07  5.74E*00 144.28E+00 0.00E+00 
i 标准 差 6.97E-56 — 7.41E-04 — 5.12E-11 1.05E-02  48.75E*00 0.00E+00 T 慰 准 差 ”5.59E-24 1.15E-05  9.27E-07  6.07E*00 41.17E+00 0.00E+00 
村 间 /s 1.2396 0.3464 0.3548 0.4984 0.7040 0.3343 时 间 /s 1.2424 0.3682 0.4162 0.5489 0.5994 0.3809 
最 优 值 0.00E+00 0.00E+00 3.66E-62  5.36E-28 2.76E-11 0.00E--00 最 优 值 3.09E-04  4.79E-04  3.10E-04  3.87E-04 3.91E-04 3.53E-04 
均值 ”3.60E-61  0.00E-00 1.33E-39 3.84E-18 6.19E-09 0.00E--00 均值  3.25E-04  2.44E-02  1.30E-03 4.7E-03 8.20E-04 6.30E-03 
i 标准 差 1.97E-60 0.00E+00 6.72E-39  1.76E-17 2.15E-08 0.00E+00 I WEZ  1.17E-05  3.39E-02  2.90E-02  8.00E-03 3.12E-04 4.20E-03 
时 间 /s 1.2610 0.3404 0.3560 0.4897 0.4967 0.3311 T/s 0.5188 0.2805 0.3443 0.2398 0.5147 0.3802 
最 优 值 0.00E+00 0.00E+00 0.00E+00 0.00E-00 3.22E-14 0.00E+00 最 优 值 -1.00E+00 -1.00E+00 -1.00E+00 -1.00E+00 1.00E+00 -1.00E+00 
均值 0.00E+00  0.00E-00  0.00E-00  3.29E-04 6.31E-10 0.00E--00 均值 -1.00E+00 -1.00E-00 -9.70E-01 -9.69E-01 -9.93E-01 -1.00E+00 
了 标准 差 0.00E+00 0.00E+00 0.00E+00 1.80E-03 9.66E-10 0.00E+00 m 示 准 差 ”0.00E+00 0.00E+00 —3.24E-02  3.19E-02 1.95E-02 5.19E-06 
时 间 /s — 1.4244 0.5261 0.5798 0.6189 0.8608 0.5910 寸 间作。 0.4247 0.1957 0.2982 0.1746 0.4395 0.5145 
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3.4 改进 蝴蝶 优化 算法 在 高 维 测 试 函 数 上 的 性 能 分 析 F9~F12，SGLBOA 的 求解 结果 均 为 理论 最 优 值 0， 而 对 于 函 

为 了 验证 SGLBOA 在 高 维 问 题 上 的 有 效 性 ,将 SGLBOA — 数 F8 mz. 同样 是 收敛 到 了 和 常见 的 最 优 值 8.88E-16。 同时 ， 
与 BOA 的 其 他 改进 算法 作对 比 。 为 保证 算法 对 比 的 时 效 性 ， 对 于 函数 F2、F8、F9~F11，SGLBOA 是 5 种 算法 中 寻 优 时 
本 文选 取 新 近 出 现 的 在 求解 函数 最 优 值 方面 性 能 较 好 的 蝴蝶 间 最 短 的 , 这 表明 SGLBOA 收敛 速度 较 快 。 对 于 函数 FA- F6 
改进 算法 作为 对 比 算法 , 即 CFSBOAP9I, MSBOA1, MSBOA2 以 及 Fl12，SGLBOA 虽然 寻 优 时 间 上 不 敌 其 他 4 种 算法 ， 但 
和 SIBOAI31。 其 中 MSBOA1 和 MSBOA2 在 对 应 的 文献 [21] 和 均值 和 方差 是 5 种 算法 中 最 小 的 ， 这 表明 SGLBOA RAE 
文献 [22] 中 都 被 记 为 MSBOA, 本 文 为 了 便于 区 分 将 其 分 别 记 为 ” 好 的 稳定 性 。 
MSBOAI 和 MSBOA2。 在 N=30 ，Max_iter=100 ， 算 法 中 的 其 由 表 6 可 知 ， 在 D=5000 时 ， 对 于 函数 F1、F2、F4~F7、 
余 参 数 与 对 应 参考 文献 中 设置 一 致 的 条 件 下 。 将 6 种 算法 分 别  F9~F12，SGLBOA 的 求解 结果 均 为 理论 最 优 值 0， 而 对 于 函 
在 D=1000、D=3000 和 DP=5000 时 独立 运行 30 次 ， 得 到 的 最 优 数 F8， 依 旧 是 收敛 到 了 常见 的 最 优 值 8.88E-16。 此 外 ， 对 于 
值 、 平 均值 、 标 准 差 、 平 均 耗 时 的 实验 数据 如 表 4~6 所 示 。 函数 F2, F8~F11, 虽然 5 种 算法 的 求解 结果 均 为 理论 最 优 值 
表 4 可知， 在 D=1000 时 ， 对 于 函数 FI-F7, F9-F12, 0, 但 SGLBOA 寻 优 时 间 最 短 。 而 对 于 函数 FI1、F4、F6、F12， 
SGLBOA 的 求解 结果 均 为 理论 最 优 值 0, 这 表明 SGLBOA 在 SGLBOA 的 寻 优 时 间 虽 然 不 是 5 种 算法 中 最 短 的 ， 但 


H 


HT 


收敛 精度 上 远 远 高 于 其 他 对 比 算法 。 而 对 于 函数 F8， 虽然 没 ” SGLBOA 的 寻 优 精度 和 稳定 性 却 是 5 种 算法 中 最 好 的 。 

有 收敛 到 理论 最 优 值 0, 但 也 收敛 到 了 常见 的 最 优 值 8.88E-16。 此 外 ， 由 表 4~6 可 知 ， 对 于 大 多 数 函 数 ， 改 进 算法 
此 外 ， 对 于 大 多 数 函数 ， 相 比 于 其 他 改进 算法 ，SGLBOA 寻 优 — CFSBOAP?, MSBOA1, MSBOA2, SIBOAP?!p]  SGLBOA 
时 间 最 短 。 因 此 SGLBOA 在 收敛 速度 上 具有 一 定 的 优势 。 相 比 于 BOA， 平 均 耗 时 都 有 所 增加 ， 但 SGLBOA 的 平均 耗 


为 了 直观 地 验证 算法 的 收敛 性 能 ， 图 5 给 出 了 6 种 算法 时 增加 量 最 小 , 而 且 SGLBOA 在 寻 优 精 度 、 收 敛 速度 、 稳 定 

在 D=1000 时 部 分 测试 函数 的 收敛 图 。 由 该 图 可 知 ,对 于 单 峰 ”性 以 及 可 扩展 性 方面 均 有 大 幅度 提升 。 
Vo HBURFI.FA-F6 以 及 多 峰 函 数 F9 和 F10, SGLBOA 在 20 代 6 给 出 了 SGLBOA 在 不 同 维 数 Dp 时 部 分 测试 函数 的 
一 左右 已 达到 收敛 ,即使 对 于 较为 复杂 的 多 峰 函 数 F11 和 F12， 收敛 图 ,通过 该 图 可 以 直观 地 看 到 随 着 测试 函数 维度 的 升 高 ， 
SGLBOA 在 40 代 左 右 也 已 收敛 ， 然 而 对 于 这 些 函 数 ， 其 他 SGLBOA 依然 保持 着 较 好 的 寻 优 性 能 。 这 表明 SGLBOA A 
对 比 算法 在 100 代 时 也 远 未 达到 收敛 ， 另外， 对 于 相同 的 欠 ” 有 良好 的 可 扩展 性 。 


代 次 数 ， 其 他 对 比 算法 的 寻 优 精度 也 都 低 于 SGLBOA。 这 表 综合 以 上 分 析 可 知 , SGLBOA 在 高 维 问题 上 的 求解 结果 

Hj SGLBOA 在 收敛 速度 和 寻 优 精度 方面 都 明显 优 于 其 余 5 整体 优 于 其 他 5 种 对 比 算法 ， 从 而 表明 SGLBOA 在 高 维 函 

种 对 比 算法 。 数 优化 问题 上 具有 较 好 的 寻 优 精度 、 收 敛 速度 、 稳 定性 和 可 
表 5 可知， 在 D=3000 时 ， 对 于 函数 F1、F2、F4~F7、 扩展 性 。 


表 4 不 同 改进 蝴蝶 优化 算法 在 D=1000 的 对 比 实验 结果 


Tab.4 Experimental results and comparison of different improved butterfly optimization algorithms with D=1000 


函数 算法 BOA CFSBOA MSBOAI MSBOA2 SIBOA SGLBOA | 函数 算法 BOA CFSBOA MSBOA1 MSBOA2 SIBOA SGLBOA 
最 优 值 2.51E-05 5.50E-83 3.93E-93 1.10E-20 4.582E-53 0.00E+00 最 优 值 1.14E-08 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 
均值 3.45E-05 3.10E-70 5.14E-17 7.72E-11 7.03E-49 0.00E+00 均值 7.01E-04 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E--00 
标准 差 4.65E-06 1.58E-69 2.82E-16 3.69E-10 2.32E-48 0.00E+00 标准 差 2.50E-03 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E--00 
= 时 间 2.4868 3.0016 3.4015 39.5388 2.6354 3.6061 时 间 /s 1.0566 1.6029 1.7189 23.0358 12591 2.0233 
PH  -— 3.75E-96 -- 1.64E-60 4.10E-04 0.00E+00 最 优 值 1.80E-03 8.88E-16 8.88E-16 8.88E-16 8.88E-16 8.88E-16 
均值 -- 2.28E-84 -- 1.14E-30 ~ 0.00E+00 均值 3.40E-03 8.88E-16 8.88E-16 9.22E-13 8.88E-16 8.88E-16 

= 标准 差  -- 6.25E-84 -- 6.20E-30 -- 0.00E--00 i 标准 差 7.49E-04 0.00E+00 0.00E--00 8.46E-13 0.00E--00 0.00E--00 
时 间 /s — -- 35.5629 313.3699 15.5259 5.1754 时 间 /s 10.7737 7.1906 13.1204 66.9593 7.7307 4.6400 
最 优 值 2.51E-19 1.05E-73 9.49E-118 2.44E+06 1.71E-58 0.00E+00 最 优 值 2.05E-05 2.28E-84 2.44E-104 4.38E-28 1.60E-61 0.00E+00 
均值 2.13E-08 4.11E-28 2.40E-21 3.57E+06 1.73E-40 1.19E+17 均值 2.32E-05 6.43E-74 1.30E-27 1.72E-06 1.01E-52 0.00E+00 

à 标准 差 4.41E-08 2.25E-27 1.32E-20 1.99E+06 9.47E-40 3.61E+17 d 标准 差 2.13E-06 3.47E-73 7.03E-27 9.43E-06 5.48E-52 0.00E--00 
时 间 /s 1.371341 1.5301 1.6819 22.5229 1.3746 1.9299 时 间 /s 17.9143 21.3937 19.1475 303.5922 17.5327 7.2015 
最 优 值 2.31E-05 1.05E-83 1.42E-98 3.56E-22 4.04E-56 0.00E--00 最 优 值 5.78E-06 0.00E--00 0.00E--00 0.00E--00 0.00E+00 0.00E--00 
均值 3.36E-05 6.29E-70 5.87E-26 3.63E-14 5.81E-50 0.00E--00 均值 2.76E-05 0.00E-00 0.00E--00 5.30E-11 0.00E--00 0.00E--00 

di 标准 差 4.34E-06 3.44E-69 3.22E-25 1.38E-13 2.59E-49 0.00E--00 dd 标准 差 1.46E-05 0.00E+00 0.00E--00 1.03E-10 0.00E--00 0.00E--00 
时 间 /s 8.0532 11.8962 8.5681 149.2748 82131 8.9430 时 间 /s 5.4525 9.5266 12.5484 107.3682 4.3755 2.3056 
最 优 值 1.93E-05 1.01E-88 3.61E-92 1.93E-51 1.00E-63 0.00E--00 最 优 值 2.05E-02 1.23E-22 1.95E-22 3.22E-13 3.01E-16 0.00E--00 
均值 2.44E-05 3.66E-78 7.88E-24 3.23E-30 4.78E-58 0.00E+00 均值 2.47E-02 7.40E-21 1.87E-07 2.86E-11 3.33E-15 0.00E+00 

i 标准 差 2.39E-06 1.57E-77 4.32E-23 1.21E-29 2.33E-57 0.00E--00 i 标准 差 1.80E-03 1.13E-20 9.84E-07 3.04E-11 2.97E-15 0.00E--00 
时 间 /s 1.0250 1.4102 1.4445 222540 1.1633 1.9457 时 间 /s 22.7010 32.6799 24.0581 399.7773 22.7013 13.5322 
最 优 值 3.20E-07 1.92E-98 2.77E-123 3.86E-80 4.02E-72 0.00E--00 最 优 值 1.80E-03 1.89E-43 1.93E-50 5.12E-22 1.00E-30 0.00E+00 
均值 1.43E-06 1.72E-85 1.63E-44 3.32E-61 7.89E-64 0.00E+00 均值 2.30E-03 4.04E-38 3.83E-11 1.92E-16 1.58E-27 0.00E+00 

m 标准 差 5.62E-07 5.94E-85 6.86E-44 1.64E-60 3.83E-63 0.00E--00 F12 标准 差 2.45E-04 1.13E-37 2.10E-10 5.54E-16 5.11E-27 0.00E+00 
时 间 /s 0.8434 1.1294 1.3750 18.2334 0.9880 1.7284 时 间 /s 2.1540 3.0340 3.6104 41.3533 2.3069 2.6784 
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表 5 不 同 改进 蝴蝶 优化 算法 在 D-3000 的 对 比 实验 结果 


Tab.5 Experimental results and comparison of different improved butterfly optimization algorithms with D=3000 


函数 ”算法 BOA CFSBOA MSBOA!1  MSBOA2 SIBOA SGLBOA | 函数 ”算法 BOA CFSBOA MSBOA!1  MSBOA2 SIBOA SGLBOA 


最 优 值 2.82E-05 2.64E-80 3.62E-97 9.58E-23 1.129E-55 0.00E+00 最 优 值 6.51E-10 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 
均值 3.43E-05 8.78E-72 1.00E-26 6.36E-14  1.75E-50 0.00E+00 均值 ”6.84E-04 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 
标准 差 3.57E-06 4.23E-71 4.13E-26 3.26E-13 8.16E-50  0.00E-00 | F7 标准 差 2.60E-03 0.00E+00  0.00E-00 0.00E+00 0.00E+00 0.00E+00 
时 间 /s ^ 2.6383 3.4132 3.7957 151.3863 4.7307 5.0902 时 间 /s ^ 2.5788 5.8771 7.3489 167.2194 3.7848 8.5168 
最 优 值 4.54E-94 5 3.06E-60 1.53E-05  0.00E--00 最 优 值 7.77E-04  8.88E-16 8.88E-16 8.88E-16 8.88E-16 8.88E-16 
均值 4.98E-82 = 1.22E-30 = 0.00E+00 均值 9.58E-04  8.88E-16 8.88E-16 3.4757e-12 — 8.88E-16 — 8.88E-16 
F2 ”标准 差 1.97E-81 6.66E-30 - 0.00E-00 | F8 标准 差 8.45E-05 0.00E+00 0.00E+00  2.2317e-12  0.00E-00 0.00E+00 
时 间 /s E 100.0776 EE 3265.1555 71.9326 27.0344 时 间 /s 23.7051 15.9489 32.6699 383.3505 20.8926 ^ 10.9917 
最 优 值 3.24E-05  1.42E-62 3.55E+08 4.24E-08 1.46E-34  4.66E-14 最 优 值 1.68E-05 4.58E-87 5.49E-93 1.54E-37 4.77E-65 0.00E+00 
均值 ”5.40E+21 3.97E-29 6.97E+18 6.98E+04 2.11E-27 5.83E+22 均值 2.10E-05 3.54E-73  1.78E-30 9.96E-16 2.15E-59 0.00E+00 


F3 标准 差 2.51E+21 2.17E-28 2.34E+19 5.13E+04 6.63E-27 1.74E+23 | F9 标准 差 2.01E-06 1.35E-72  9.74E-30 4.47E-15 8.04E-59 — 0.00E--00 


时 间 /s 5.717 5.0131 7.7128 187.3159 5.7349 7.6199 时 间 /s 51.1329 54.9721 45.666 2074.5216 42.5611 25.3653 
最 优 值 ”3.04E-05  1.35E-79  6.18E-88 5.18E-22 3.12E-55  0.00E--00 最 优 值 1.09E-06 0.00E+00  0.00E-00 — 0.00E-00  0.00E-00  0.00E--00 
均值 3.57E-05  4.55E-67  9.90E-19 8.37E-12 1.09E-49  0.00E--00 均值 6.62E-06 0.00E-00  0.00E-00  2.4253e-13  0.00E-00 0.00E+00 


FA 标准 差 3.23E-06  2.49E-66 — 5.42E-18 4.57E-11 3.87E-49  0.00E-00 || FIO 标准 差 5.62E-06  0.00E-00  0.00E-00  1.3284e-12  0.00E-00 0.00E+00 


时 间 /s 38.4395 60.7089 45.0255 2030.0856 40.2137 41.6311 时 间 /s 23.6051 22.8525 31.2439 399.7311 16.2036 8.673 
最 优 值 2.05E-05  5.36E-90 1.24E-89 9.86E-45 1.21E-63 ”0.00E+00 最 优 值 2.03E-02 8.23E-24  4.52E-19 6.11E-14 3.73E-16 — 0.00E-00 
均值 2.58E-05  3.84E-74 5.01E-27 5.06E-28 2.26E-57 0.00E+00 均值 2.47E-02 7.39E-21 2.93E-09 5.31E-11 3.21E-15 0.00E+00 


F5 标准 差 2.66E-06 2.09E-73 2.75E-26 2.35E-27 9.33E-57 ”0.00E+00 | F11 标准 差 2.06E-03 1.15E-20 1.32E-08 9.47E-11 4.63E-15 0.00E+00 


时 间 /s 2.7723 5.1965 5.8508 192.7266 3.1643 7.8005 时 间 /s 72.4257 102.9056 74.4172 3543.2641 64.372657 40.5479 
最 优 值 9.48E-09  8.56E-99 1.08E-108 4.36E-77 8.99E-73 0.00E+00 最 优 值 1.93E-03  8.84E-46 9.38E-43 1.65E-22 7.57E-31 0.00E+00 
均值 5.34E-07 1.68E-83  9.83E-48 1.55E-64 7.58E-68 0.00E+00 均值 2.28E-03 2.15E-38  244E-14 4.31E-17 1.57E-28 0.00E+00 


F6 标准 差 2.46E-07 9.17E-83  3.75E-47 ”8.02E-64 2.21E-67  0.00E-00 | F12 标准 差 1.64E-04 5.31E-38 1.34E-13 1.05E-16 4.51E-28 0.00E+00 


时 间 /s ^ 3.1369 4.1799 5.4877 142.1456 3.9304 7.0109 时 间 /s ^ 5.4693 8.0107 9.4143 244.6403 5.7093 7.7365 


表 6 不 同 改进 蝴蝶 优化 算法 在 D=5000 的 对 比 实验 结果 


Tab.6 Experimental results and comparison of different improved butterfly optimization algorithms with D=5000 


函数 算法 BOA CFSBOA MSBOAI MSBOA2 SIBOA SGLBOA| 函 数 算法 BOA CFSBOA MSBOAI MSBOA2 SIBOA  SGLBOA 
最 优 值 2.77E-05 2.15E-81 6.12E-90 8.47E-24 1.11E-54 0.00E+00 最 优 值 2.85E-11 0.00E+00 0.00E+00 0.00E+00  0.00E--00 0.00E+00 
- 均值 3.43E-05 3.37E-68 9.38E-21 6.10E-13 1.94E-51 0.00E--00 - 均值 ”3.21E-08 0.00E--00 0.00E-00 0.00E--00 0.00E+00 0.00E--00 
标准 差 3.31E-06 1.84E-67 5.14E-20 3.28E-12 5.71E-51 0.00E--00 标准 差 4.09E-08 0.00E--00 0.00E--00 0.00E--00 0.00E+00 0.00E--00 
时 间 /s 7.7651 8.4849 8.8395 391.4922 7.3067 10.1121 时 间 /s 5.9967 9.2575 10.2146 324.8185 6.8590 . 10.2753 
最 优 值 -- 5.29E-93 -- 8.03E-58 6.42E-05 0.00E+00 最 优 值 5.53E-04 8.88E-16 8.88E-16 8.88E-16  8.88E-16  8.88E-16 
b 均值 -- 8.17E-80 B 1.96E-25 -  0.00E+00 -" ^M 6.46E-04 8.88E-16 8.88E-16 8.88E-16  8.88E-16  8.88E-16 
标准 差 。 -- 4.36E-79 = 1.07E-24 --  0.00E+00 示 准 差 5.59E-05 0.00E+00 0.00E+00 3.87E-12 0.00E+00 0.00E+00 
时 间 /s -- 99.4272 ET 8296.3997 83.2912 27.489 时 间 /s 32.9034 24.0397 44.565 718.5489 25.1074 14.9316 
最 优 值 2.59E-05 1.20E-61 6.81E-08 3.81E*04 4.16E-31 3.86E+18 最 优 值 1.00E-05 3.03E-87 4.41E-96 8.69E-30  7.95E-64 0.00E+00 
m 均值 1,94E+22 3.04E-26 145E-19 1.37E+05 7.58E-25 1.28E425 io 均值 2.38E-05 4.47E-74 9.18E-40 2.74E-14  3.53E-57 0.00E+00 
标准 差 3.78E+22 1.16E-25 4.51E+19 7.81E+04 1.69E-24 1.99E425 示 准 差 2.68E-06 2.34E-73 4.34E-39 1.25E-13  1.90E-56 0.00E+00 
时 间 /s 9.6822 9.1985 11.8947 400.2919 9.2538 9.4833 时 间 /s 64.5646 75.5385 63.5889 4959.0585 60.8643 24.8883 
最 优 值 2.87E-05 2.12E-82 2.09E-83 6.55E-21 6.39E-54 0.00E--00 最 优 值 3.72E-06 0.00E--00 0.00E--00 0.00E--00  0.00E--00 0.00E--00 
- 均值 3.60E-05 1.45E-69 2.43E-17 5.65E-11 3.44E-49 0.00E--00 im ^[É 1.93E-05 0.00E+00 0.00E-00 2.47E-11 0.00E+00 0.00E--00 
标准 差 4.22E-06 7.91E-69 1.33E-16 3.09E-10 1.26E-48 0.00E--00 7.47E-06 0.00E+00 0.00E--00 8.46E-11 0.00E+00 0.00E--00 
时 间 /s 38.8434 55.9658 40.7198 3044.8077 38.9418 42.0427 34.003 33.3223 42.9037 1312.8956 25.3992 13.2231 
最 优 值 1.91E-05 1.63E-88 2.64E-84 4.90E-58 1.88E-62 0.00E--00 1.33E-02 7.47E-23 1.71E-19 6.61E-14  1.93E-16 0.00E+00 
均值 2.69E-05 4.85E-76 7.29E-16 1.85E-28 2.09E-56 0.00E--00 1.46E-02 9.04E-21 3.04E-08 1.18E-11  2.79E-15 0.00E+00 
id 标准 差 3.21E-06 1.74E-75 3.99E-15 8.33E-28 9.86E-56 0.00E--00 5.93E-04 1.32bE-20 1.64E-07 1.67E-11 3.02E-15 0.00E+00 
时 间 /s 7.32330 9.0121 | 92273 4122919 7.7815 10.4181 时 间 /s 102.4434 145.7607 105.7969 9219.1753 188.959827 57.0347 
最 优 值 1.86E-08 7.17E-97 2.21E-106 1.31E-90 2.82E-73 0.00E--00 最 优 值 2.09E-03 2.49E-45 1.50E-40 2.40E-23  5.57E-31 0.00E+00 
a 均值 7.05E-07 7.09E-84 2.90E-36 6.36E-65 1.80E-66 0.00E+00 its 均值 2.55E-03 3.66E-38 3.35E-1l] 749E-17  5.98E-28 0.00E+00 
标准 差 3.24E-07 3.86E-83 1.49E-35 2.15E-64 9.32E-66 0.00E--00 标准 差 2.10E-04 1.39E-37 1.78E-10 1.85E-16  2.12E-27 0.00E+00 


时 间 /s 6.1523 7.3095 8.8976 | 280.3555 6.3400 | 9.2357 时 间 /s 13.6142 18.798 23.1732 884.6742 14.8423 17.9933 
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Fig.6 Convergence diagrams of SGLBOA on some test functions in 


3.5 Wilcoxon 秩 和 检验 


对 于 改进 算法 ， 通 常 需要 利用 统计 检验 来 验证 ， 相 对 于 


现 有 算法 而 言 ， 改 进 算 法 是 否 具 


有 显著 性 。 为 判断 SGLBOA 
HJ 


的 寻 优 结果 是 否 在 统计 上 与 


岂 算 法 的 最 佳 结果 显著 不 同 ， 


different dimensions 


性 显著 ;否则 就 是 接受 HO 假设 ， 说 明 两 种 算法 在 整体 上 是 


(rae 
给 出 1 


相同 的 。 表 7 
BOAB!, CFSBOAPU, 


在 14 个 基准 测试 函数 下 ，SGLBOA 对 
MSBOA1、MSBOA2 以 及 SIBOAD3] 的 


秩 和 检验 ? 值 ,从 表 7 中 可 以 观察 到 大 多 数 的 值 是 小 于 5%, 


Wilcoxon 秩 和 检验 在 5% 的 显著 性 水 平 下 进行 也 即 当 P 值 小 姑 此 总 体 上 说 明 SGLBOA 与 其 他 5 种 算法 之 间 差 异性 显著 ， 
于 5% 时 ,，“R" 为 “+”， 即 拒绝 HO 假设 ,说 明 两 种 算法 的 差异 ”SGLBOA 明显 优 于 其 他 5 种 对 比 算法 。 
表 7 Wilcoxon 秩 和 检验 结果 
Tab.7 Wilcoxon rank sum test results 

BOA CFSBOA | MSBOA] | MSBOA2 SIBOA BOA CFSBOA | MSBOAI | MSBOA2 SIBOA 

p R p R p R p Rp R p R p Rp Rp R p R 
Fl  12118e-12 + 12118e-12 + 12118e-12 + 1.2118e-12 + 12118e-12 -| F8  12118e-12 + NA = 03337  - 86979e-07 + NA - 
F2  12118e-12 + 12118e-12 + 12118e-12 + 12118e-12 + 12118e-12 -| F9  12118e-12 + 1.2118e-12 + 1.2118e-12 + 12118e-12 + 12118e-12 + 
F3  12118e-12 + 12118e-12 + 12118e-12 + 12118e-12 + 12118e-12 +| FI0  12118e-12 + NA = NA = 00056 - NA = 
F4  12118e-12 + 1.2118e-12 + 1.2118e-12 + 1.2118e-12 + 12118e-12 +| FII  12118e-12 + 1.2118e-12 + 1.2118e-12 + 1.2118e-12 + 12118e-12 + 
F5  12118e-12 + 12118e-12 + 1.2118e-12 + 1.2118e-12 + 1.2118e-12 +| F12  12118e-12 + 1.2118e-12 + 1.2118e-12 + 1.2118e-12 + 12118e-12 + 
F6  12118e-12 + 12118e-12 + 1.2118e-12 + 1.2118e-12 + 1.2118e-12 +| F13 7.0430e-07 + 1.8567e-09 + 1.6351e-05 + 8.9934e-11 + 2.4386e-09 + 
F7 1.2118e-12 + NA = NA = NA i NA =|| F14 42770e12 + NA - NA B NA - NA - 


注 :“R” 为 显著 性 判断 结果 ， 


4 ARA 


本 文 在 蝴蝶 优化 算法 的 基础 上 提出 拉丁 超 立 方 抽样 的 自 


适应 高 斯 小 孔 成 像 蝴 蝶 优 化 算法 。 通 过 将 概率 统计 中 的 拉丁 


“+”“-”“=” 分 别 表 示 MSBOA 的 性 能 优 于 、 劣 于 和 相当 对 比 算 法 ; 


“NA” 表 示 不 适 


， 无 法 进行 显著 判断 


超 立 方 抽样 策略 应 / 


] 于 算法 的 种 群 初始 化 ， 提 高 种 群 的 多 样 


性 ， 


以 提升 算法 的 全 局 搜索 能 力 ; 


适应 最 优 引导 策略 的 引 


入 能 够 自动 调节 不 同 进化 时 
局 部 搜索 能 力 ， 从 而 提升 算法 的 寻 优 精度 ; 引入 高 


的 全 局 和 


期 的 搜索 范围 ， 较 好 地 平衡 算法 


录用 定稿 R A, F: 
斯 小 孔 成 像 策 略 步 提升 算法 的 寻 优 精度 并 加 快 算法 的 


KIGE. 14 个 测试 函数 的 仿真 实验 结果 表明 SGLBOA 算 
法 具有 较 好 的 寻 优 精度 、 收 敛 速度 、 稳 定性 以 及 可 扩展 性 。 
同时 由 于 SGLBOA 算法 在 高 维 优化 问题 上 具有 较 好 的 性 能 ， 
未 来 可 进一步 研究 将 算法 应 用 于 大 规模 优化 问题 中 。 
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